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Abstract 

Electrochemical  and  electrical  double  layer  characteristics  for  the  phase  boundary  nanoporous  carbon  1 1 .0  M  triethylmethylammonium 
tetrafluoroborate  (TEMA)  solution  in  acetonitrile  (AN),  7-butyrolactone  (GBL),  acetone  (DMK)  and  propylene  carbonate  (PC)  have  been 
studied  using  the  cyclic  voltammetry  (CV)  and  the  electrochemical  impedance  spectroscopy  (EIS)  methods.  The  constant  current  and 
constant  power  charge/discharge  characteristics  of  capacitor  modules  constructed  and  filled  with  different  electrolytes  were  tested  over  the 
temperature  range  from  —30  to  +60  °C.  Tests  for  “lifetime”  of  the  capacitors  were  performed  at  elevated  temperature.  Energy  density  versus 
power  density,  i.e.  so-called  Ragone  plots  were  constructed  to  characterize  the  performance  of  the  capacitor  modules.  It  was  established  that 
the  energy  density,  power  density,  discharge  time,  etc.  of  the  electrical  double  layer  capacitors  (EDLC)  depend  on  the  pore  size  distribution 
and  conductivity  of  the  nanoporous  electrode  material.  It  was  shown  that  the  energy  output  of  capacitor  depends  on  the  viscosity  and  molar 
conductivity  of  a  solvent.  The  power  output  of  the  capacitor  increased  in  the  order  of  solvent  PC  <  DMK  <  GBL  <  AN  in  agreement 
with  the  dependence  of  the  Ragone  plot  parameters  on  temperature. 
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1.  Theoretical  background 

Electrical  double  layer  capacitors  (EDLC)  (supercapaci¬ 
tors)  are  the  intermediate  systems  between  electrochemical 
batteries  that  can  store  high  energy  density  associated  with 
low  power  density  values,  and  dielectric  capacitors,  which 
can  deliver  very  high  power  during  few  milliseconds  [1].  Su¬ 
percapacitors  can  generate  high  specific  power  during  few 
seconds  or  more,  which  leads  to  the  specific  energy  densi¬ 
ties  ranging  from  0.5  to  lOWhkg-1.  There  are  three  dif¬ 
ferent  types  of  supercapacitors  based  on  the  different  elec¬ 
trode  materials  and  working  principles:  (1)  carbon/carbon, 
(2)  metal  oxide  and  (3)  electronically  conducting  polymers. 
This  work  is  focused  on  carbon/carbon  supercapacitors  [1]. 

The  important  advantages  of  carbon/carbon  EDLC  are 
their  very  good  reversibility  and  the  comparatively  low  tem¬ 
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perature  coefficient.  As  the  electrostatic  interactions  are  sig¬ 
nificantly  less  detrimental  to  electrodes  and  solution  stabil¬ 
ity  than  the  usual  electrochemical  redox  reactions,  used  for 
the  generation  of  electricity  in  the  fuel  cells  and  various  bat¬ 
teries,  EDLC  can  be  recharged-discharged  up  to  106  times 
[1].  It  is  well  known  that  the  performance  characteristics  of 
an  EDLC,  e.g.  in  terms  of  the  relations  between  achieved 
power  densities  and  corresponding  energy  densities  depend 
on  the  equivalent  series  resistance  (ESR)  and  the  internal 
distribution  of  electrode  resistance  (IER)  in  the  pore  matrix 
of  the  electrodes  [1-17]. 

The  capacitance  of  EDLC  depends  mainly  on  the  carbon 
material  used  for  preparation  of  the  electrodes.  Theoreti¬ 
cally,  the  higher  is  the  surface  area  of  the  activated  carbon, 
the  higher  specific  capacitance  should  be  expected.  How¬ 
ever,  the  practical  situation  is  more  complicated  and  usually 
the  capacitance  measured  does  not  have  the  linear  relation¬ 
ship  with  the  specific  surface  area  of  the  electrode  material. 
The  main  reason  for  this  phenomenon  is  that  the  nanopores 
with  small  diameter  may  not  be  accessible  to  the  electrolyte 
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solution  simply  because  the  electrolyte  ions,  especially  big 
organic  ions  and  ions  with  the  solvation  cell,  are  too  big 
to  enter  into  the  nanopores.  Thus,  the  surface  area  of  these 
non-accessible  nanopores  will  not  contribute  to  the  total  dou¬ 
ble  layer  capacitance  of  the  electrode  material. 

The  most  frequently  employed  methods  for  testing  of  the 
EDL  capacitors  are  (1)  the  cyclic  voltammetry  (CV)  for  gen¬ 
erating  the  charge-discharge  curves;  (2)  galvanostatic  meth¬ 
ods  such  as,  e.g.  current  step,  and  constant  current  meth¬ 
ods  (CC)  at  I  =  constant;  (3)  constant  power  method  at 
P  =  constant  (CP);  (4)  the  electrochemical  impedance  spec¬ 
troscopy  (EIS). 

1.1.  Cyclic  voltammetry 


In  the  case  of  CV  applied  at  a  cell  potential  sweep-rate 
v  =  ±dU/dt  (Vs-1),  the  current  response  of  a  capacitor 
with  the  capacitance  C  is  given  as 


dU 

dt 


Cv 


Quite  often,  for  carbon  double  layer  capacitors,  Faradaic 
pseudocapacitance  devices  and  hybrid  devices  C  is  a  func¬ 
tion  of  cell  potential  U  during  charge  or  discharge  cycle,  i.e. 
C(U)  =  f(U).  Then 

dU 

i  =  C(U ) — =  C(U)v  (2) 

dt 


Thus,  from  Eqs.  (1)  and  (2)  C(U)  is  directly  obtainable  as 
i(U)/v ,  where  i(U)  is  a  potential  and  potential  sweep  rate 
dependent  current  density.  It  should  be  noted  that  by  CV  at 
moderate  scan-rates  a  differential  capacitance  C  =  dq/dU  is 
directly  generated  rather  than  the  integral  capacitance  C  = 
Aq/AU.  The  integral  between  potentials  U\  and  U2,  in  CV 
gives  directly  the  charge  A q,  passed  between  these  two  po¬ 
tentials: 


C(U)dU 


1.2.  Electrochemical  impedance  spectroscopy 


Electrochemical  impedance  spectroscopy  is  based  on  a 
modulation  function  U(t)  =  Vo  sin  cot,  where  co  =  Inf  is 
the  angular  frequency  and  /  the  ac  frequency  (0.01  Hz  < 
/  <  100  kHz).  The  current  response  of  the  electrical  double 
layer  (approximated  as  a  capacitor  and  resistor  in  series)  is 
given  as 

dU(t) 

i(t)  ~  C - =  C coo  Vo  sin  (cot  +  cp)  (4) 

dt 

where  (p  is  a  phase-shift  that  equals  to  90°  for  the  ideal 
capacitor  and  Vo  the  alternating  voltage  signal  amplitude. 
For  electrochemical  systems  usually  cp  <  90°,  depending  on 
the  presence  of  any  effective  or  equivalent  series  resistance 
or  on  any  slowness  of  polarization  of  the  dielectric  material 
of  the  capacitor.  In  electrochemical  double  layer  capacitors 


the  polarization  slowness  arises  from  the  dipole  orientation 
of  the  solvent  molecules  and/or  from  the  ion  adsorption,  but 
it  is  observable  only  at  higher  frequencies  [1-3]. 

The  response  function  in  EIS  for  an  ideal  capacitor  can 
be  simply  represented  by  imaginary  part  of  the  impedance 
Z",  given  as 


In  real  EDLC  systems,  Z"  is  normally  coupled  with  an  ohmic 
resistance  R  that  is  either  an  actual  or  an  ESR  and  corre¬ 
sponding  impedance  is  simply  written  as 

Z'  =  R  (6) 


Thus,  for  a  real  EDL  capacitor  (series  circuit) 

Z  =  Zf  +  Z"  =  R+^—  (7) 

J  coC 

from  which  R  and  C  can  be  determined  by  frequency-response 
analysis  in  well-known  ways  [1-3]. 

For  the  real  electrochemical  double  layer  capacitor,  both 
R  and  C  are  quantities  which  are  distributed  over  a  range 
of  the  effective  values  due  to  the  electrode  porosity,  i.e. 
nanoporous  structure  of  the  electrode,  electrolyte,  EDLC 
separator  parameter,  design,  etc.  and  the  impedance  Z  as  a 
function  of  ac  frequency  is  a  more  complicated  parameter 
than  that  in  Eq.  (7).  According  to  de  Levie  [8,10,11]  and 
Salitra  et  al.  [5]  the  real  system  has  a  spectrum  of  “RC” 
time-constants  and  a  corresponding  “power  spectrum”.  The 
complex  impedance  plane-diagram,  i.e.  the  —  Z"  versus  Z! 
dependence  over  a  range  of  frequencies  consists  mainly  of 
two  principal  regions:  one  at  45°  line  at  high  frequencies, 
extrapolating  to  an  intercept  on  the  Z! -axis  equal  to  the  exte¬ 
rior  solution  +  constant  resistances;  the  other  —  a  transition 
at  a  “knee”  in  the  curve  to  a  vertical  line  (ideal  capacitor) 
for  which  —  Zl'  is  simply  j  /coC  at  sufficiently  low  frequen¬ 
cies.  The  extrapolation  of  this  line  to  the  Zl- axis  gives  the 
total  distributed  and  external  resistance  values.  The  so-called 
pore  resistance  Rp0re  is  defined  as  the  real  component  of  the 
impedance  when  co  ->  0. 


1.3.  Constant  current  study 


In  the  current  step  technique  a  constant  current  1  is  applied 
on  the  cell  and  the  potential  response  is  measured.  For  a 
capacitive  system  with  an  ohmic  resistance  the  potential  drop 
A  U  at  very  short  times  r  can  be  expressed  as 

A  U  =  IR  (8) 


In  the  constant  current  study  (CC)  on  charging  curve  (U, 
t- curve  at  I  =  constant)  A q  is  directly  determinable  from 
the  record  of  U(t)  versus  t ,  i.e.  U(t )  versus  q.  The  capaci¬ 
tance  of  the  electrode  material,  as  well  as  the  cell  can  be 
evaluated  from  the  slope  of  the  discharge  curves  according 
to  the  formula  below: 
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When  for  simplicity  the  charging  curve  is  approximated  by 
a  linear  function  in  the  region  d  U,  the  medium  integral  ca¬ 
pacitance  values  can  be  obtained  as 


At 

I - 

A  U 


(10) 


From  the  CC  curves  the  internal  resistance  can  be  evaluated 
from  IR-drop  after  the  current  is  changed.  Regularly  the 
response  time  10  ms  is  used,  where  dU\  is  the  value  of  the 
voltage  during  10  ms  and  I  the  current  used: 


(11) 


1.4.  Constant  power  study 


For  double  layer  capacitors  the  peak  power  capability  can 
be  estimated  from  a  constant  current  discharge  from  the 
rated  voltage  Uw  to  Uw/2  according  to  the  Ref.  [1].  Thus, 
the  average  voltage  for  the  discharge  is  (3/4 )UW.  At  these 
conditions  the  average  efficiency  EF  for  the  discharge  is 


EF  =  1  —  / 


R 


(3/4)  Uw 


(12) 


The  current  I  for  the  discharge  at  an  efficiency  EF  is  defined 
as 

(3/4)(l-EF  )UW 


I  = 


R 


(13) 


and  the  corresponding  power  is  expressed  as 
9(1 -EF  )UW 


P  = 


16  R 


(14) 


The  matched  power  density  is  expressed  as 


^matched  — 


a 


w 


4  R 


(15) 


The  power  density  (PD)  of  the  electrical  double  layer  ca¬ 
pacitor  can  be  calculated  using  Eq.  (16): 

PD  =  -  (16) 

m 


where  m  is  the  weight  of  the  capacitor. 

The  power  density  of  electrical  double  layer  capacitor  for 
a  specified  efficiency  is  given  by 


PD  (W  kg- 1 )  = 


(17) 


where  Ki  =  C/m  [Ah(kgcell)-1]  (when  the  cell  capaci- 
tance  is  expressed  in  ampere  hours,  Ah)  and  K\  =  CR.  The 
energy  density  ED  is  expressed  as 


ED  (Wh  kg-1)  = 


k2uI 

7200 


(18) 


It  should  be  noted  that  the  discussion  given  above  is  valid  for 
EDLC,  regularly  the  constant  power  method  has  been  used 
and  the  capacitor  has  been  cycled  between  maximum  cell 
potential  Uw  and  Uw/2  and  the  power  densities  are  between 
0.1  and  10  Wg"1. 


2.  Experimental 

A  nanoporous  carbon  powder  was  made  from  titanium 
carbide  (A.  C.  Starck,  CA  grade,  1-3  jxm)  by  a  chlorinating 
method  [18].  TiC  with  an  average  particle  size  of  1 .3-3.0  jxm 
was  loaded  into  the  silica  rotary  kiln  rector  and  let  react 
with  a  flow  of  chlorine  gas  (99.999%  pure)  for  4h  in  a  tube 
furnace  at  950  °C.  The  by-product  (TiCLj.)  was  led  away 
by  the  stream  of  the  excess  CI2.  After  reaction  the  reactor 
was  flushed  with  the  argon  (99.9995%  pure)  at  1000  °C  for 
0.5  h  to  remove  the  excess  CI2  and  residues  of  a  gaseous 
by-products  from  carbon.  During  heating  and  cooling  the  re¬ 
actor  was  flushed  with  a  slow  stream  of  argon.  Thereafter  the 
resulting  carbon  powder  was  moved  into  silica  stationary  bed 
reactor  and  treated  with  hydrogen  at  800  °C  for  2.5  h.  Dur¬ 
ing  heating  and  cooling  the  reactor  was  flushed  with  a  slow 
stream  of  helium  (99.9995%  pure).  Carbonaceous  polariz¬ 
able  electrodes  from  these  nanoporous  carbon  powders  were 
prepared  as  follows.  The  mixture  of  94%  (wt.)  nanoporous 
carbon  and  6%  (wt.)  polytetra  fluoroethylene  PTFE  binder 
(Aldrich,  60%  suspension  in  water)  was  rolled  stepwise  into 
the  carbon  film  with  a  final  thickness  of  100  db  5  jxm.  After 
drying,  the  raw  electrode  sheets  were  plated  from  one  side 
with  a  thin  aluminum  layer  (4  ±  1  |xm)  using  the  plasma 
activated  physical  vapor  deposition  method. 

It  should  be  noted  that  the  carbon  powders  with  different 
pore  size  distribution,  shown  in  Fig.  1,  were  used  to  produce 
a  positively  and  negatively  charged  electrodes.  Methods  for 
modification  of  the  pore  structure  to  achieve  a  better  match 
between  the  pore  size  and  the  electrolyte  ion  size  are  partic¬ 
ularly  discussed  in  Ref.  [4,15-18].  Some  characteristics  of 
these  carbon  powders  such  as  specific  surface  area,  microp¬ 
ore  volume  and  area,  obtained  from  low  temperature  nitro¬ 
gen  sorption  data  using  ASAP  2000  (Micromeritics,  Inc.) 
analyzer,  are  presented  in  Table  1 .  Single-point  BET  surface 
area  has  been  calculated  at  relative  pressure  P/Po  =  0.2 
and  the  total  pore  volume  at  near  to  a  saturation  pressure, 
P/Po  =  0.99.  The  volume  and  surface  area  of  micropores 
were  derived  from  the  so-called  Fplot  dependence  of  ni¬ 
trogen  sorption  isotherm  using  the  Harkins-Jura  statistical 


Fig.  1.  Pore  size  distribution  of  carbon  powders  1  and  2  according  to  the 
density  functional  theory. 
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Table  1 


Pore  characteristics  of  nanoporous  carbon  powders  used  for  the  preparation 
of  electrodes 


Parameter 

Carbon  no.  1 

Carbon  no.  2 

BET  surface  area,  Sa  (m2g_1) 

1449 

1922 

Micropore  surface  area,  Sm  (m2g_1) 

1518 

1749 

Micropore  volume,  Vm  (cm3  g~ 1 ) 

0.63 

0.83 

Total  pore  volume,  Vtot  (cm3  g-1) 

0.80 

1.09 

Median  pore  diameter3,  Dhk  (A) 

6.8 

9.4 

a  Median  pore  diameter  according  to  Horvath-Kawazoe  method. 


thicknesses.  The  pore  size  distribution  was  calculated  ac¬ 
cording  to  the  density  functional  theory  (DFT). 

The  electric  double  layer  capacitors  (EDLC)  were  as¬ 
sembled  so  that  the  positivly  charged  electrodes  were  made 
from  the  carbon  no.  1  with  a  micropore  predominantly  be- 

o 

low  10  A  and  the  negatively  charged  electrodes  from  the 
carbon  no.  2  with  a  significant  fraction  of  micropores  be- 

o 

tween  10  and  30  A.  The  electrodes  were  attached  to  current 
collector:  Al  foil,  and  interleaved  with  an  ion-permeable 
separator  paper  from  Codashi  Nippon.  The  electrode  pairs 
from  positively  and  negatively  charged  electrodes  with  a 
visible  surface  area  of  64.5  cm2  were  connected  in  parallel. 
The  electrode  pack  prepared  was  placed  in  a  sealed  box, 
kept  at  100  °C  under  vacuum  for  3  days  to  remove  all  gases 
absorbed  and  then  impregnated  with  electrolyte.  Four  EDLC 
cells  were  filled  with  different  electrolytes  comprising 
1.0  M  triethylmethylammonium  tetrafluoroborate  (TEMA) 
in  anhydrous  acetonitrile  (AN,  Riedel-de  Haen,  E^O  < 
0.003%),  anhydrous  Selectipur®  propylene  carbonate  (PC, 
Merck,  H2O  <  0.003%),  Selectipur®  7-butyrolactone 
(GBL,  Merck,  H2O  <  0.003%)  and  Selectipur®  dimethyl 
ketone  (DMK,  Alfa,  H20  <  0.03%).  TEMA  salt  (Stella) 
was  dried  under  vacuum  at  170  °C  prior  to  preparing  the 
electrolyte  solution.  The  triethylmethylammonium  tetraflu¬ 
oroborate  has  been  selected  out  because  the  solubility  of 
this  electrolyte  is  noticeable  higher  at  lower  temperatures  in 
comparison  with  tetraethylammonium  tetrafluoroborate  and 
therefore  the  more  higher  conductivities  for  electrolyte  at 
lower  temperatures  can  be  achieved.  Selected  physical  char¬ 
acteristics  of  the  solvents  used  in  the  EDLCs  are  presented 
in  Table  2  [19]. 

2.7.  Electrochemical  study 

The  EDLCs  were  preconditioned  by  continuous  cycling 
between  2.3  and  1.15  V  with  the  current  /  =  1  A  during  3 


days  prior  performing  the  further  electrochemical  studies. 
It  was  observed  that  approximately  1000  cycles  (that  took 
about  24  h)  were  needed  to  stabilize  the  inner  resistance  and 
capacitance  of  AN  based  EDLC.  The  EDLCs  based  on  PC, 
GBL  and  DMK  electrolytes  demonstrated  stable  characteris¬ 
tics  of  polarization  after  2000  cycles  (approximately  2  days). 
After  preconditioning  the  pressure  release  valve  was  opened 
to  remove  the  gases  formed  during  continuous  potential  cy¬ 
cling. 

2.2.  The  cyclic  voltammetry  (j,  E)  curves 

The  cyclic  voltammetry  curves  for  EDLCs,  based  on 
1.0  M  TEMA  in  AN,  GBL,  DMK  and  PC  electrolyte  are 
presented  in  Fig.  2.  The  EDLCs  with  the  AN,  GBL  and  PC 
electrolytes  are  ideally  polarizable  in  the  potential  region 
from  0  to  3  V.  The  experimental  data  revealed  that  for  all 
the  EDLCs  the  j ,  ^-curves  are  independent  of  the  number 
of  current  cycles,  n,  if  n  >  5. 

At  higher  potential  scan  rates,  the  well  expressed  distor¬ 
tion  effects  observed  in  the  /,  ^-curves  for  GBL,  DMK  and 
PC  are  caused  by  the  higher  internal  resistance  of  the  EDLCs 
compared  with  the  acetonitrile  based  electrolyte  EDLCs,  as 
well  as  by  the  onset  of  so-called  “electrolyte  starvation”  ef¬ 
fect,  discussed  by  Conway  and  Pell  [1,16]  and  associated 
with  the  withdrawal  of  the  electrolyte  ions  from  the  pore 
bulk  electrolyte  due  to  the  adsorption  of  ions  on  the  double 
layer  interfaces  when  they  become  charged. 

In  Fig.  3a  and  b,  it  is  seen  that  the  current  responses  ex¬ 
pectedly  increase  with  higher  scan  rate  values.  The  voltam- 
mograms  become  distorted  from  the  ideal  rectangular  form 
characteristic  for  an  ideal  capacitor  demonstrating  constant 
capacitance  [1].  The  corresponding  capacitance  behavior  for 
EDLCs  are  shown  in  Fig.  3c,  which  indicates  that  the  capac¬ 
itor  based  on  1.0  M  TEMA  in  acetonitrile  (Fig.  3c)  behaves 
as  almost  ideal  EDLC  already  at  scan  rate  of  50mVs_1. 
The  GBL  based  capacitor  behaves  quite  similarly  to  AN 
system,  however  the  deviation  from  the  rectangular  shape  at 
sweep  rate  of  50  mV  s_1  is  obvious.  The  capacitor  modules 
consisting  acetone  and  propylene  carbonate  electrolyte  so¬ 
lutions  (Fig.  3d)  behave  as  ideal  capacitors  at  scan  rates  be¬ 
low  10  mV  s  .  Furthermore,  slightly  increased  capacitance 
in  the  module  with  acetone  at  cell  potential  values  exceed¬ 
ing  2.0  V  are  indication  of  the  more  narrow  electrochemical 
window  (i.e.  of  the  more  narrow  region  of  ideal  polarizabil¬ 
ity)  and  probably  indicates  the  beginning  of  the  Faradaic 
reactions  in  the  EDLC.  The  increased  capacitance  at  higher 


Table  2 

The  physical  parameters  of  the  solvents  [19] 


Solvent 

Melting  point  (°C) 

Boiling  point  (°C) 

Viscosity  (Pas 

Dielectric  constant  ( s ) 

Density  (kgm  3) 

AN 

-43.8 

81.6 

0.369 

36.64 

0.786 

GBL 

-43.3 

204.0 

1.72 

39.0 

1.128 

PC 

-48.8 

241.7 

2.513 

66.14 

1.205 

DMK 

-94.8 

56.0 

0.306 

21.01 

0.790 
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Fig.  2.  The  cyclic  voltammetry  curves  at  potential  scan  rate  of  10  mV  s  1  for  the  EDLCs  in  1.0  M  TEMA  in  AN,  GBL,  DMK  and  PC. 


voltage  is  decrease  of  the  effective  size  of  electrolyte  ions 
at  increased  voltage.  This  effect  was  recently  observed  and 
discussed  for  the  skeleton  carbon  electrodes  [2,3]. 

2.3.  Galvano static  cycling  study  of  the  cells 

The  EDLCs  were  tested  at  constant  current  charge/discharge 
regimes  in  voltage  range  from  1.15  to  2.3  V.  The  current 
was  varied  from  0.5  to  8.0  A.  The  discharge  capacitance 
was  calculated  from  the  data  of  third  cycle  according  to 


Eq.  (10).  Measurements  were  repeated  at  various  fixed  tem¬ 
peratures  in  the  region  —  30  °C  <  T  <  +60  °C  to  study  the 
stability  of  EDLCs  parameters.  The  internal  resistance  value 
was  calculated  from  the  IR — drop  (at  fixed  d t  =  10  ms) 
according  to  Eq.  (11). 

The  shortest  discharge  cycle  observed  for  EDLCs  based 
on  DMK,  presented  in  Fig.  4,  corresponds  to  the  lowest  ca¬ 
pacitance  value  of  EDLC  modules  studied.  The  similar  ef¬ 
fect  is  also  observed  at  low  current  values  (Fig.  5).  It  should 
be  noted  that  the  capacitance  decreases  almost  proportion- 


Fig.  3.  Current  density  vs.  capacitor  potential  curves  for  AN  (a),  and  PC  (b)  electrolyte  corresponding  EDLC  systems  at  potential  scan  rates  of  50,  10 
and  5  mV  s-1  noted  in  figure.  Respective  C,  E-curves  (c,  d)  calculated  from  j,  E- curves  at  potential  scan  rates  of  50,  10  and  5  mV  s-1,  noted  in  figure. 
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Fig.  4.  Constant  current  charge/discharge  cycle  at  I  =  0.5  A  and 
T  =  298  K  for  the  EDLCs  filled  with  different  electrolytes,  noted  in  figure. 


Fig.  5.  Constant  current  discharge  profiles  measured  for  EDLC  at 
T  =  298  K  filled  with  1.0  M  TEMA  in  AN  at  constant  current  densities, 
noted  in  figure. 

ally  with  increasing  the  current  with  exception  of  PC  based 
EDLC  (Fig.  6).  The  dependence  of  C  versus  I  increases  with 
decreasing  temperature  and  is  also  dependent  on  the  elec¬ 
trolyte,  particularly  on  viscosity  of  the  solvent.  The  specific 
capacitance  and  resistance  data  are  presented  in  Table  3. 

The  capacitance  values  presented  in  Table  3  decrease  from 
AN  to  PC  electrolyte  at  ambient  temperatures.  The  tempera¬ 
ture  dependence  is  most  significant  for  PC-based  EDLC  that 


Fig.  6.  Discharge  capacitance  vs.  current  density  plots  at  T  =  298  K  for 
the  EDLCs  based  on  different  nonaqueous  electrolytes,  noted  in  figure. 


Table  3 


The  specific  capacitance  and  resistance  values  for  the  EDLC  cells 


Solvent 

Capacitance  (F) 

Resistance  (£2  cm2) 

-30  °C 

+20  °C 

+60  °C 

— 30°C 

+20  °C 

+60  °C 

AN 

22.4 

22.7 

22.9 

1.0 

0.6 

0.5 

GBL 

19.8 

21.6 

22.7 

3.0 

1.2 

1.0 

DMK 

17.6 

21.0 

21.7 

4.5 

1.2 

1.3 

PC 

6.1 

19.7 

21.5 

7.9 

1.8 

1.2 

The  capacitance  for  the  cells  is  calculated  according  to  Eq.  (10)  and  the 
resistance  per  visible  electrode  surface  area  according  to  Eq.  (11). 


is  also  well  supported  by  the  highest  viscosity  among  the 
solvents  studied.  Consequently  the  logical  result  was  the  al¬ 
most  stable  capacitance  and  resistance  over  the  temperature 
range  from  —30  to  +60  °C  for  the  EDLC  comprising  ace¬ 
tonitrile  as  the  low-viscosity  solvent.  The  specific  resistance 
value  (in  £2  cm2)  is  influenced  by  the  choice  of  EDLC  ma¬ 
terials  (i.e.  carbon,  separator,  current  outlet,  etc.)  and  design 
of  the  EDLC  cell.  As  a  rule,  the  smaller  is  the  capacitance 
value,  the  lower  are  the  specific  series  resistance  obtained. 
In  present  work,  all  cells  had  the  same  design,  therefore  the 
resistance  values  should  be  mainly  influenced  by  the  proper¬ 
ties  of  the  electrolyte  and  carbon.  However,  according  to  the 
systematic  analysis  of  experimental  results,  the  resistance 
values  correlate  well  with  the  electrolyte  molar  conductivity 
values. 

The  cycling  efficiency,  so-called  round  trip  efficiency 
(RTE)  can  be  calculated  as  the  ratio  of  capacitances  mea¬ 
sured  during  discharging  and  charging  of  EDLC,  i.e.  RTE  = 
C(discharge)/C(charge),  where  the  discharge  and  charge 
capacitance  values  are  calculated  according  to  Eq.  (10). 

Among  capacitors  studied,  the  round  trip  efficiency  is  the 
highest  in  acetonitrile-based  EDLC  and  the  RTE  value  de¬ 
creases  in  the  order  of  solvents  AN  >  GBL  >  DMK  >  PC 
(Table  4).  RTE  is  also  influenced  by  temperature:  the  lower  is 
the  temperature,  the  lower  is  RTE.  EDLC  with  acetone-based 
electrolyte  has  a  RTE  maximum  at  ambient  temperature  that 
could  serve  as  indication  of  increased  irreversibility  of  the 
processes  in  DMK  based  EDLC  at  elevated  temperature. 
RTE  value  generally  decreases  when  the  current  density  is 
increased,  but  with  exception  of  DMK.  Although  the  latter 
phenomenon  of  acetone-based  electrolyte  is  not  clear  yet, 
it  could  be  caused  by  the  worse  cycleability  of  ketone  type 
solvents  in  the  wide  region  of  electrode  polarization. 


Table  4 

Round  trip  efficiency  (RTE)  for  EDLC  cells  with  different  electrolytes 
measured  at  variable  temperature  and  current  density  (/  =  0.5  and  2.0  A) 


Solvent 

RTE  (%)  at  /  =  0.5  A 

RTE  (%)  at  I  =  2.0  A 

-30  °C 

20  °C 

+60  °C 

-30  °C 

20  °C 

+60  °C 

AN 

95.9 

97.6 

97.9 

87.4 

93.2 

94.9 

GBL 

86.1 

93.9 

96.3 

70.9 

86.4 

90.7 

DMK 

83.1 

92.2 

80.6 

68.9 

85.6 

83.9 

PC 

73.5 

90.2 

93.1 

63.6 

76.5 

82.1 
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Table  5 

Capacitance  in  farads  of  the  EDLC  calculated  using  data  of  different 
testing  methods 


Solvent  Measurements  procedure  for  capacitance  (F) 


CV 

( v  =  lOmV  s-1) 

CC 

(7  =  0.5  A) 

EIS  at  /  =  10  mHz 
(DC  =  2.3  V, 

AC  =  5  mV) 

AN 

22.9 

22.7 

25.6 

GBF 

23.1 

22.8 

26.0 

PC 

21.9 

21.6 

22.8 

DMK 

20.4 

19.5 

23.9 

According  to  the  data  in  Table  5,  the  capacitance  val¬ 
ues  of  the  EDLC  cells  depend  of  the  measurement  tech¬ 
nique  used.  The  DC-capacitance  values  obtained  from  CC 
or  CV  measurements  are  almost  similar,  but  substantially 
different  from  those  calculated  from  the  impedance  data. 
In  all  cases  the  impedance  spectrum  yields  the  higher  se¬ 
ries  capacitance  Cs  values  for  charged  EDLCs.  The  excep¬ 
tional  case  is  the  PC-based  EDLC,  which  yields  the  lower 
Cs  values — probably  because  the  used  frequency  10  mHz 
is  still  too  high  to  establish  the  equilibrium  Cs  capacitance 
value  in  the  high-viscosity  PC  electrolyte. 

2.4.  Electrochemical  impedance  study 

From  the  Z",  7!  plots  (Fig.  7)  it  can  be  found  that  all  the 
capacitors  behave  as  typical  capacitors  over  the  frequency 
range  from  100  Hz  to  10  mHz.  The  curve  in  figure  can  be 
divided  into  three  main  sections.  At  high  ac  penetrability  (at 
low  frequency  /  <  0.1  Hz)  the  carbon  material  behaves  like 
porous  material  and  the  ac  signal  detects  the  large  amount  of 
the  pore  volume.  This  region  of  Z",  7!  plot  is  called  a  planar 
section  and  the  phase-angle  approaches  asymptotically  to 
-80  •••  -90°. 

The  second  section  lies  in  the  region  10  Hz  <  /  <  0.1  Hz, 
where  ac  signal  detects  the  less  porous  material  and  the 


Fig.  7.  Nyquist  plots  for  the  EDFC  based  on  different  solvents,  noted  in 
figure. 


curve  slope  on  Cole-Cole  or  complex  plane  plots  (Fig.  7)  is 
almost  45°.  The  characteristics  of  plateau  with  the  45°  slope 
depend  on  many  factors  including  carbon  properties,  elec¬ 
trode  thickness,  separator  and  electrolyte  characteristics, 
temperature,  etc.  EDLC  cells  under  this  study  differed  from 
each  other  only  by  the  choice  of  the  solvent.  Therefore,  the 
differences  observed  in  the  “45° -slope”  region  character¬ 
istics  are  mainly  caused  by  the  different  properties  of  the 
electrolyte  solution  and  may  particularly  be  contributed  with 
the  electrolyte  (solvent)  viscosity.  In  the  high-frequency 
region  of  ac  the  porous  electrode  behaves  like  a  flat  surface 
and  therefore  the  capacitance  obtained  is  by  many  orders  of 
magnitude  lower  compared  to  the  equilibrium  capacitance  at 
co  — >  0.  This  region  of  Z",  Z'-plots  (usually  at  /  =  100  Hz) 
is  used  to  determine  the  series  resistance  of  the  double-layer 
capacitor,  Z'^_^  =  Rmt.-  The  resistance  component  ob¬ 
served  at  7"  =  0  in  Fig.  7  satisfactorily  matches  the 
electrolyte  conductivity  and  resistance  values  presented  in 
Tables  3  and  6. 

Based  on  the  data  in  Table  6  it  was  concluded  that  the 
specific  capacitance  (Fern-2)  is  weakly  influenced  by  the 
choice  of  the  solvent  in  electrolyte  composition  as  long  as 
the  same  concentration  of  electrolyte  salt  is  considered.  In 
the  case  of  propylene  carbonate  the  somewhat  lower  Cs  value 
is  probably  due  to  the  high  viscosity  of  PC  that  does  not 
allow  reaching  the  limiting  capacitance  at  relatively  high 
frequency  (/  =  10  mHz).  The  data  in  Tables  1  and  6  well 
confirm  the  influence  of  electrolyte  molar  conductivity  on 
the  so-called  pore  resistance  and  electrolyte  resistance  val¬ 
ues.  The  volumetric  capacitance  of  the  carbon  electrodes 
has  almost  similar  value  in  all  EDLC  cells  tested  in  this 
study.  Our  investigations  reveal  the  clear  advantage  of  ace¬ 
tonitrile  based  capacitor,  particularly  due  to  the  low  total 
impedance. 

Relationships  between  series  resistance  and  series  capac¬ 
itance  versus  frequency  are  presented  in  Fig.  8  and  well  sup¬ 
port  the  data  demonstrated  in  Table  6.  In  Fig.  8,  it  is  also 
seen  that  all  test  cells  have  a  plateau  on  the  Cs -curve  at  very 
low  ac  frequency  that  is  characteristic  of  the  EDL  formation 
in  fine  micropores. 

Fig.  9  shows  the  phase  angle  versus  frequency  relation¬ 
ship  for  all  electrolyte  compositions  under  study,  however 
the  ideal  behavior  of  the  capacitor  according  to  the  phase 
angle  is  not  reached  in  any  electrolyte  studied.  Perhaps  the 
closest  to  the  ideal  one  is  the  cell  based  on  AN  with  the  phase 
angle  value  of  8  ~  —90°,  which  indicates  the  nearly  pure 
capacitive  behavior  of  EDLC  at  low  frequencies.  The  cells 
filled  with  GBL  and  DMK  electrolytes  have  slightly  lower 
absolute  phase  angle  values  compared  to  the  AN  electrolyte, 
which  may  be  caused  by  the  specific  interactions  between 
solvent  bearing  carbonyl  group  and  nanoporous  carbon  [1]. 
The  same  dependence  seems  to  be  true  for  the  PC,  how¬ 
ever  the  much  lower  phase  angle  value  measured  at  constant 
frequency  indicates  the  existence  of  the  additional  limiting 
factors,  one  of  which  could  be  the  low  molar  conductivity 
of  PC  electrolyte. 
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Table  6 

Electrochemical  parameters  of  the  EDLC  cells  calculated  from  electrochemical  impedance  spectroscopy  data  at  2.3  V  (DC) 


Solvent 

Cs  (Fern  2) 

Re i  (£2  cm2) 

^pore  (^2  Cm  ) 

Rs  at  /  =  100  Hz 
(Q  cm2) 

Average  capacitance  of  carbon  electrodes 

(Fg-1)  (F  cm-3) 

AN 

0.40 

0.29 

1.20 

0.31 

106 

71 

GBL 

0.40 

0.65 

2.79 

0.75 

106 

69 

DMK 

0.37 

0.70 

3.91 

0.77 

102 

70 

PC 

0.35 

1.37 

4.94 

1.67 

94 

65 

(a)  ft  Hz  (b)  f/  Hz 

Fig.  8.  Series  resistance  and  series  capacitance  of  the  cells  based  on  different  solvents,  noted  in  figure. 


Table  7 


Energy  and  power  performance  of  the  EDLC  cells,  calculated  from  galvanostatic  cycling  test  data  at  the  discharge  potentials  up  2.3-1.15  V 


Solvent 

Usable  energy 
(Wh  dm-3) 

Rs  (Q  cm2) 

Power  densitya 
(90%  eff)  (W cm-3) 

Real  max.  powerb 
(W  cm-3) 

Matched  powerc 
(W  cm-3) 

AN 

5.5 

0.63 

6.3 

15.5 

59.6 

GBL 

5.5 

1.18 

3.6 

6.5 

31.8 

DMK 

5.4 

1.21 

3.6 

6.5 

32.1 

PC 

4.5 

1.81 

2.2 

2.6 

19.5 

a  Power  density  at  90%  efficiency  according  to  Eqs.  (14)  and  (16). 
b  Real  maximum  power  obtained  from  Ragone  plot  (Fig.  10)  when  E  ->  0. 
c  Matched  power  according  to  Eq.  (15). 


2.5.  Constant  power  (CP)  study 

Fig.  10  presents  the  Ragone  plot,  i.e.  the  energy  and  power 
relationship  of  the  unpacked  capacitors  that  is  derived  from 
CP  curves.  Here,  the  unpacked  capacitor  is  defined  as  the 


ft  Hz 


EDLC  cell  without  case.  More  precisely,  it  includes  the  vol¬ 
ume  and  weight  of  terminals,  carbon  electrodes,  separator 
and  electrolyte  solution. 

The  volumetric  energy  densities  at  moderately  low  power 
density  were  rather  similar  in  all  capacitors.  This  is  in  a  good 
accordance  with  the  capacitance  values  obtained  from  con¬ 
stant  current  (CC)  study.  The  behaviour  of  capacitor  filled 
with  acetonitrile  retains  almost  constant  energy  in  the  wide 


Fig.  9.  The  plot  of  phase  angle  vs.  frequency  for  EDLC  cells  with  different  Fig.  10.  Ragone  plot  for  the  unpacked  EDLC  cells  with  different  electrolyte 
solvents,  noted  in  figure.  solutions,  measured  at  +20  °C. 
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range  of  studied  power.  The  capacitors  filled  with  acetone 
and  y-butyrolactone  electrolytes  behave  rather  similarly:  the 
significant  energy  decrease  at  higher  power  values  was  ob¬ 
served.  Similar  behavior  of  these  two  electrolytes  was  also 
seen  in  other  figures  (Figs.  6-9).  The  behaviour  of  capaci¬ 
tor  filled  with  propylene  carbonate  was  significantly  differ¬ 
ent  from  other  studied  cells — it  had  noticeable  energy /power 
dependence  already  at  rather  low  power  values. 

The  maximum  power  values  (at  E  =  0)  calculated  from 
the  Ragone  plots  and  presented  in  Table  7  well  agree  with  the 
internal  resistance  values  of  capacitors  (cf.  data  in  Tables  3 
and  6).  Power  density  values  in  Table  7  calculated  according 
to  Eqs.  (14)  and  (16)  match  with  the  experimental  values 
measured  at  RTE  90%. 

3.  Conclusions 

The  electrical  double  layer  capacitor  (EDLC)  cells  com¬ 
prising  nanoporous  carbon  electrodes  and  different  elec¬ 
trolyte  solutions — 1.0M  triethylmethylammonium  tetraflu- 
oroborate  in  acetonitrile  (AN),  y-butyrolactone  (y-BL), 
acetone  (DMK)  and  propylene  carbonate  (PC) — were  thor¬ 
oughly  studied  using  several  electrochemical  testing  proce¬ 
dures  and  techniques.  It  was  found  that  most  of  important 
parameters  such  as  the  series  capacitance,  series  resistance 
and  phase  angle  noticeably  depend  on  the  electrolyte  com¬ 
position  as  well  as  on  the  cell  potential  applied.  It  was 
established  that  the  so-called  pore  resistance,  i.e.  internal 
distribution  of  the  electrode  resistance  rises  in  the  following 
sequence  of  solvent:  AN  <  GBL  <  DMK  <  PC.  It  was 
shown  that  the  value  of  characteristic  frequency  depends 
on  the  electrolyte  solution  and  on  the  electrode  potential, 
i.e.  on  the  nature  and  structure  of  solvated  ions  adsorbed  at 
the  surface  of  the  nanoporous  carbon  electrode.  It  should 
be  noted  that  the  results  of  the  galvanostatic  tests  and 
estimated  Ragone  plots  show  the  clear  advantage  of  the 
acetonitrile-based  electrolytes  in  carbon/carbon  EDLCs. 
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